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Abstract 

We arrange the loopwise perturbation theory for the critical vis- 
cosity exponent Xr^, which happens to be very smah, as a power series 
in Xr^ itself and argue that the effect of loops beyond two is negligible. 
We claim that the critical viscosity exponent should be very closely 
approximated by Xrj = + ^) — 0.0685. 
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Critical exponents, amplitude ratios and scaling functions were issues of 
considerable importance three decades ago. Sophisticated calculations and 
experiments were carried out which clearly established the correctness of the 
various theoretical models (Landau-Ginzburg equations for statics and the 
various models of dynamics [1-4] introduced by Hohenberg and Halperin, 
Kawasaki, Ferrell etc). Basically, the exponents could be classified into two 
types: i) large exponents i.e. exponents of 0(1) and ii) small exponents i.e. 
exponents of 0(0. 1) or even smaller. It is the small exponents where the 
most crucial confrontation between theory and experiment can occur. That 
is why even after three decades the small exponents remain an interesting 
issue. In static critical phenomena [5] the small exponents are associated 
with the critical correlation function at the transition point (77, the anoma- 
lous dimension exponent) and specific heat {a, the specific heat exponent, 
the specific at constant volume for the liquid-gas transition and the specific 
heat at constant pressure for the superfiuid transition of He'^), while in the 
critical dynamics the small exponent is associated with the shear viscosity. 
Accurate determination [6] of a for the superfiuid transition and comparison 
with very detailed calculation [7] confirm the theoretical expectation. For the 
shear viscosity exponent z^, the recent measurements [8] in the space shuttle 
have yielded an accurate value, namely 2;^ = 0.0690 ±0.0006. The theoretical 
self consistent 2-loop calculation in D=3 of Hao yields Zj^ — 0.066 ± 0.002 , 
amazingly close to the experimental value. This raises the immediate ques- 
tion: what happened to the higher loops? The one loop answer is 20 % away 
from the experimental answer, the two loop calculation produces the 20 % 
enhancement almost entirely and so what happens to the infinite number 
of loops that have been left out. This is the question that we address in 
the paper and provide the insight into why the higher loops happen to be 
unimportant. 

In a liquid-gas system near the critical point or a binary liquid mixture 
near the critical mixing point, the order parameter is the density (concen- 
tration) difference and relaxes when disturbed from equilibrium according to 
the Langevin equation 

^M. = _rk\e + K')m + N{k) (1) 

where (f){k) is the Fourier transformation of the D-dimensional field (f){xi, ..xd)- 
In the relaxation rate the factor fc^ indicates that (i>{x) is conserved. F is 
the Onsagar coefficient and the diffusion constant is D=~; where x is the 
susceptibility. Near the critical point the susceptibility is x = (^^ + n"^)'^ 
with K = the inverse correlation length which diverges near T = Tc as 
^ oc \T — Tc\~''. The term N is a stochastic forcing that comes from the short 
wavelength modes. Fluctuation dissipation holds and the correlation of N is 
related the usual way to the dissipation. 

In a fluid the density (concentrations) fluctuations will be affected by the 
velocity fluctuations and the effect of the velocity is to advect the concentra- 
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tion field, so that 



^ + tka E v^imk - P) = -rk\k' + K^)m + N{k) (2) 
p 

The fact that the velocity fluctuations affect the concentration means that 
we need to know the velocity fluctuations. The equation of motion (for small 
fiuctuation) is Navier-Stokes equation 

dva{k) ,2 



dt 



-rjk'vM + Kik) (3) 



Note Va and Na are solenoidal. However Eqs.(2) and (3) do not conserve the 
local free energy density J2k[{k'^ + n^)4>{k)(j){—k) + v{k)v{—k)] when the dis- 
sipation terms are omitted and consequently Eq.(3) needs to be argumented 
as 

dva{k) 



dt 



+ iY.P%T^p{^)cl){p)(i){k -p) = -r^k^M + Kik) (4) 



Where Tap{k) — Sap — the projection operator. The effect of the non- 
linear terms in Eqs. (4) and (5) is to renormalize the Onsager co-efficient F 
and the shear viscosity rj. Dropping the non-linear terms, we get the zeroth 
order solution 

0(0) t)= I e-^^"^^"+^"^^'-''^N{t')dt' (5) 

and 

v^^\k,t)= / e-'^^"^'-''^N^{t')dt' (6) 



The first order solution is easily seen to be 

0^^) t) = -j e-^^^(^^+«^)(*-*')ifc. E iP^ {k - P, t')dt' (7) 

p 

and 

v'^^\k, t)^-i e-''*^'(*-*'Vp;30(°)(P, - t')dt' (8) 

The fields being stochastic in nature, the effect of the non-linear terms in 
Eq.(2) and Eq.(4) are to be understood as averaged over the noise terms 
and it is easy to see that the non-linear terms in Eq.(2) yields a term of 
the form -k'^{k'^ + k^) ^T^^\k,t- t')<p'-°\-k,t')dt' and those in Eq.(4) give 
—k'^ J r]^^\k,t — t')v^\—k,t')dt' , when we split the quadratically non-linear 
term as one field at zeroth order and the other at first order. For Eq.(2) 
this imphes writing the non-hnear term as ika[< '}lpV^a\p)^^^\k — p) > -|- < 
J2pV^Hp)4>^°Kk — p) >] and similarly for Eq.(4). This is exactly equivalent 
to a one loop result. The two loop results come from all the pairings of 3, the 
three loop from the pairings of 5 and so on. The Fourier transforms r^{k, uj) 
and ri^{k,uo) of T^{k,t — t') and 'r]^{k,t — f ) are the renormalized Onsager 
co-efficient and the shear viscosity respectively. 
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The renormalized transport coefficients T^^\k,uj) and ri^^\k,uj) diverge 
at the critical point in the zero frequency, zero wavelength limit and dom- 
inates the molecular contributions. From now on we will refer to these as 
T{k,uj) and ri{k,uj). A little algebra shows that at one loop, we get the 
standard results {ri{k)k^ > r{k)k'^{k^ + k^)) 

, If d^p sin^O 

and 

r]{k,K) 



(I'p ir{jr — p'^) sin 9 



4C3 J (p2 + K2)(p/2 + ^2) ^2(^2 + ^2)r(p) + p'^p'^ + K^W)] 

(10) 

where p' = k — p. 

We now introduce the scaling behavior (long wavelength divergence at 
the critical point) at k = as 

T{k) - roA;-^+^'' (11) 

r^{k) = r)ok-''^k^ (12) 

consistent with Eqs.(9) and (10), where is the exponent that is yet un- 
known. Working at k = in Eq.(9), we find 

r/oro = y + 0(x,) (13) 

We anticipate at this stage that is very small and are going to use it as a 
small parameter in setting up our calculation. Our main observation is that 
a loopwise expansion can be cast as an expansion in powers of Xr, for the 
quantity rjoTo. We can get yet another expansion for ryoFo by using Eq.(lO) 
at K = 0. The integral has a long wavelength divergence at = and this 
leads to the evaluation of the integral as a pole in x^. This yields 



VoTo = + 0(1) (14) 



Combing Eqs(13) and (14), 

= ll^ (15) 

to the lowest order. 

We now observe that the perturbation theory for T and 77 can be expressed 
through diagrams as shown in Fig.l and 2, with a wavy line denoting the 
velocity field (propagator or correlator as the case may be) and a solid line 
the density field. 

For higher loops, self energy insertion [13-15] are not shown separately. 
They are handled by a finite frequency evaluation of a lower loop and yields 



4 



insignificant corrections. The important graphs are the vertex correction 
varieties [16] that are shown in Fig.l and 2. If one compare a one loop and 
a two loop graph, we note that compared to the one loop graph, the two 
loop graph has two time zones, one dominated by viscosity relaxation, the 
other lacking any viscosity contribution. This means an additional factor of 
(^oro)~^ everytime a loop increases. Now, in addition we note that for every 
loop the viscosity graphs diverge logarithmically if = and has a pole 
for small x^. We simply need to evaluate this pole in a manner very similar 
to the dimensional regularization scheme in field theory. Prom Fig.l, there 
emerges 

Vo^o = Ji + -^J2 + rTM^-^s + (16) 

while from Fig. 2, we get 

Vo^o = - + + . + (17) 

where and J„ are integrals of which the one loop parts, Ji and Ji. Fig. la 
and 2b are shown in Eqs.(9)and (10). The integrals corresponding to the two 
loop (Fig. lb and 2b) 



1 r d'y f 

2 J 'C^ J 



W C's p^q^^k — p — q)'^{p'^ + {k — fp}^){q^ + {k — q)'^) 

(18) 



Jo 



21 



dPp r (fq [{k — p — q)"^ ~ (k — p)'^][{k — p — q)"^ ~ {k — q) 

J C3 (J^ — jpj'iQ^ _ j^2(J^ _ ^ _ ^2 

kaTap{p){k - q)pkf,Tf,y{q){k - p)i, 
p^q^[{k — p)^ + {k — q)^ + {k — p — q)^] 

(19) 

Using Eq.(14) to substitute for t^qFo in Eqs.(16) and (17), we end up with 
TyoFo = Jl + 15Xr,J2 + (15x^)^3 + (20) 

The important fact is that 15x^72 -C 1 and this trend continues through 
higher loops. This fixes t/qFo- 

Turning now to the diagrams of Fig.2, they lead to (using Eq.(14) repeat- 
edly) 



ryoFo = — + I5/2 + x,(15)% + 



Xjj 



— [1 + 15Xr,^ + Xrf{15)'^Xrf^ + 



1 M 8 8 15/3 2 8 152/4 

-— 1 + - + x,-—^ + x,'-—^ + 21 

15a;,, TT^ h h 
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leading to the ordering in Xj,. The calculation of ^ yields | and hence to two 
loop order 

^^^ll^^^ + 3^)^°-°^^^ ^^^^ 
The reason why I2 is smaller than Ii has to do with the projection factors 
which yield zeroes in the integrand. The large number of zeroes and their 
distributions in the three loop integral leads to ^ being significantly smaller 
than The additional factor of now makes the three loop contribution 
negligible. The important point is that for a n-loop integral /„, the projection 
factor produce sufficient cancellation that IS^^^/n is always of 0(1) and that 
ensures that higher loops produce insignificant corrections when <^ 1. 

The generic form of the three loop graph of Fig. 2c involves a few different 
time ordering, all of which are shown in Fig. 3. A typical contribution 13^^ 
coming from the last two graphs of Fig.(3)is 

^2^(1) f (Pp d^q d^r [{k — p}"^ — p^] [{k — r)^ — r^] 

C3 C's C's {k — p)'^{k — q)'^ 

PaTai3{.k)rpp^T^^{p - q)p^ 



X- 



X 



{k - r)2[p3 + ^|3][g3 ^\k- g]3] 

QjT^xil - r)q\ 

[r^ + \k — f\^]{k — p — q)'^{k — q — r)' 



(23) 

The evaluation of J3 has to be in the limit of ^ 0. This allows us to 
drop 'A;' from all the the terms after a factor of k"^ has been extracted from 
the integral. We now carry out the following steps in a D-dimensional space 
for generality 

a) expand the number in powers of A;^ and keep the first term (this is 
proportional to A;^) and set k=0 everywhere else. 

b) do an angular average over the directions of k. 

In a D-dimensional space, l'^^ after some long algebra reduces to 



^3 - 4^(^ + 2) y \^{q- T)\p- q)^ 

\p^q^-{p-q)'^][q^r'^ - (g • r)^] 



X- 



(24) 



„2^2 



It is the factor {p-rf — which is qualitatively new. The two loop integral 
I2 did not have such a factor. The characteristic feature of this factor is 
that in the absence of the quite indirect additional appearance of the angle 
between p and r because of the term \q — r|^and \p — in the denominator, 
the averaging over the directions of r (or p) would make l'^^ identically zero. 



6 



In practice, this effect makes it unusually small compared to I2 or Ii, which 
do not have such a factor. If we look at the higher loops, each additional loop 
brings in a factor of this type and that is the reason behind the successive 
diminishing of each of these integrals. 

A numerical evaluation yields I^^^ ~ "(l^)^ which makes the point that 
wanted to make. The correction from the three loop graphs are down by an 
order of and this effect persists to higher orders. This is the reason why the 
two loop calculation of the viscosity exponent gives an answer surprisingly 
close to the experimental value. 

In closing we would like to mention that we have used a gaussian free 
energy in this calculation. There is a quartic part in the free energy which 
is responsible for the anomalous dimension rj. The correction coming from 
this is once again largest at the loop level when it first appears. Higher loop 
graphs involving the four point vertex give a much smaller contribution once 
again because of the frequent zeroes in the integrand at higher loops. The 
net result is that from one to two loops there is a substantial change in 
but thereafter the contribution of the higher loops are ordered by itself 
and with the integrals themselves quite small, the small value of ensures 
that the higher loop effects are small. 
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FIGURE - CAPTIONS 



Fig.l : One and two loop diagrams for the density relaxation rate. The 
solid lines stand for density fluctuation and the wavy lines for velocity 
fluctuation. 



Fig. 2 : One and two loop diagrams for the viscosity. The solid lines 
stand for density fluctuation and the wavy lines for velocity fluctuation. 



Fig. 3 : Three loop diagrams of the vertex correction variety for the 
viscosity showing all possible time orderings. Propagators appear with 
an arrow and correlators with an open circle. 
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